Bragg spectroscopy of a Bose— Einstein condensate 



Os 

a 

(N 



> 
OS 

o 



o 

os 



i 

c 

o 
o 



13 



J. Stenger, S. Inouye, A. P. Chikkatur, D.M. Stamper-Kurn, D.E. Pritchard, and W. Ketterle 
Department of Physics and Research Laboratory of Electronics, 
Massachusetts Institute of Technology, Cambridge, MA 02139 

February 1, 2008 

Properties of a Bose-Einstein condensate were studied by stimulated, two-photon Bragg scattering. 
The high momentum and energy resolution of this method allowed a spectroscopic measurement of 
the mean-field energy and of the intrinsic momentum uncertainty of the condensate. The coherence 
length of the condensate was shown to be equal to its size. Bragg spectroscopy can be used to 
determine the dynamic structure factor over a wide range of energy and momentum transfers. 



The first evidence for Bose-Einstein condensation in 
dilute gases was obtained by a sudden narrowing of the 
velocity distribution as observed for ballistically expand- 
ing clouds of atoms Q. Indeed, most textbooks describe 
Bose-Einstein condensation as condensation in momen- 
tum space However, the dominant contribution to 
the observed momentum distribution of the expanding 
condensate was the released interaction energy (mean- 
field energy) resulting in momentum distributions much 
broader than the zero-point motion of the ground state 
of the harmonic trapping potential. Since the size of a 
trapped condensate with repulsive interactions is larger 
than the trap ground state, the momentum distribution 
should be considerably narrower than in the trap ground 
state. In this paper, we show the momentum distribution 
of a trapped condensate to be Heisenberg uncertainty 
limited by its finite size. This is equivalent to showing 
that the coherence length of the condensate is equal to 
its physical size. 

Sub-recoil momentum resolution has been previously 
achieved by resolving the Doppler width of a Raman tran- 
sition between different hyperfine states || or of a two- 
photon transition to a metastable excited state Q . Here 
we use Bragg scattering where two momentum states of 
the same ground state are connected by a stimulated 
two-photon process p|. This process can be used to 
probe density fluctuations of the system and thus to mea- 
sure directly the dynamic structure factor S(q,^), which 
is the Fourier transform of the density-density correla- 
tion function and is central to the theoretical descrip- 
tion of many-body systems Q . In contrast to measuring 
S(q,f ) with inelastic neutron scattering like in superfluid 
helium M , or using inelastic light scattering ||||] , Bragg 
scattering as used here is a stimulated process which 
greatly enhances resolution and sensitivity. 

Bragg scattering of atoms from a light grating was 
first demonstrated in 1988 Jl(J and has been used to ma- 
nipulate atomic samples in atom interferometers jnj, in 
de Broglie wave frequency shifters |Pj, and also to cou- 
ple out or manipulate a Bose-Einstein condensate 
Small angle Bragg scattering, called recoil-induced res- 



onances, has been used for thermometry of laser-cooled 
atoms jl3| . In this work we establish Bragg scattering as 
a spectroscopic technique to probe properties of the con- 
densate. We refer to it as Bragg spectroscopy in analogy 
to Raman spectroscopy which involves different internal 
states. 

The absorption of N photons from one laser beam 
and stimulated emission into a second laser beam con- 
stitutes an iV-th order Bragg scattering process. The 
momentum transfer q and energy transfer hv are given 
by q = 2Nhks'm('d/2) and v = NAv, where $ is the an- 
gle between the two laser beams with wave vector k and 
frequency difference Av. 

For non-interacting atoms with initial momentum Tiki , 
the resonance is given by the Bragg condition hv — 
q 2 /2m + hkiq/m, which simply reflects energy and mo- 
mentum conservation for a free particle. The second term 
is the Doppler shift of the resonance and allows the use 
of Bragg resonances in a velocity-selective way [| ] 
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For a weakly interacting homogeneous condensate at 
density n, the dispersion relation has the Bogoliubov 
form j|] 



^Jv 2 + 2v nU/h, 



(1) 



where nU = nAirh 2 a/m is the chemical potential, with 
a and m denoting the scattering length and the mass, 
respectively, and hvo = q 2 /2m. At low energies the exci- 
tation spectrum is phonon-like, obeying hv ~ cq, where 
c is the speed of sound 0. For energies hv » nU the 
spectrum is particle-like 



v w vq + nU /h. 



(2) 



The mean-field shift nU reflects the exchange term in 
the interatomic interactions: a particle with momentum 
q experiences twice the mean field energy as a particle 
in the condensate || . We use this property to determine 
the condensate mean-field energy spectroscopically. This 
is related to, but different from the mean-field shift due 
to interactions with an electronically excited state which 
was used to identify BEC in atomic hydrogen . 
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Eq. (Q) is the excitation spectrum of a homogeneous 
condensate with initial momentum hki = 0. The inho- 
mogeneous trapping potential adds two features which 
broaden the resonance: a distribution of initial momenta 
due to the finite size of the cloud, and a distribution of 
mean-field shifts due the density variation over space. 

The momentum distribution along the x-axis is given 
by the Fourier transform of the wavefunction: f(p x ) — 
[f dx dy dz e~ w ' cX / h '$>(x,y, z)] . In the Thomas-Fermi 
approximation the wave function ty(x,y, z) in an har- 
monic trapping potential is [^S(x, y, z)\ = no(l — 
(x/xo) 2 — (y/yo) 2 — (z/zo) 2 ), where no denotes the peak 
density. The size of the wavefunction is determined by 
no and the trapping frequencies vi through the coeffi- 
cients £o,?/o, z o : == \f^ n aU /m(2m> x ) 2 (and similar 
for yo,za). The line shape I p {p x ) is proportional to the 
square of the Fourier coefficients [[L4| 



IpiPx) ~ [J2(PxX /h)/(p x x /h)' 2 



(3) 



where J 2 denotes the Bessel function of order 2. This 
curve is very similar to a Gaussian and has an rms-width 
of Ap x = y/21/8 Ti/xq. Thus, the corresponding Doppler 
broadening Av p = y/21/8 q/2-KmxQ of the Bragg reso- 
nance is inversely proportional to the condensate size Xq 
and does not depend explicitly on the number of atoms. 

The same parabolic wavefunction gives the (normal- 
ized) density distribution (15n/4no)\/l — n/n Q . The 
simplest model for the spectroscopic lineshape {v) due 
to the inhomogenous density assumes that a volume el- 
ement with local density n leads to a lineshift nil (eq. 



15h(i/ — v$) 
4n U 



'1- 



n U 



(4) 



which has a total width of noU/h, a maximum at 
2noU/3h, an average value (first moment) of An^U/lh, 
and an rms-width of Av n = y / o/lA7noU/h. In con- 
trast to the Doppler broadening due to the finite size, 
the mean-field broadening only depends on the density 
but not explicitly on the size. 

In our experiments the combined broadening mecha- 
nisms represented by eqs. ([|) and (Q) have to be con- 
sidered. While the exact calculation of the lineshape 
requires detailed knowledge of the excitation wavefunc- 
tions, the total line shift (first moment) and rms-width 
can be calculated using sum rules and Fermi's Golden 
Rule. Thus, it can be rigorously shown that the to- 
tal line shift remains 4nof7/7/i, while the rms-widths 

Av = ^/Ai/^ + is the quadrature sum of the Doppler 

and mean-field widths jlq]. 

We produced magnetically trapped, cigar-shaped 
Bose-Einstein condensates as in previous work jl(|. In 
order to study the resonance as a function of density and 



size, we prepared condensates using two different trap- 
ping frequencies and varied the number of atoms by re- 
moving a variable fraction using the rf-output coupler 
]l7t . The density of the condensate was determined from 
the expansion of the cloud in time-of-flight and the size 
from the atom number and the trapping frequencies [jl6| . 
Bragg scattering was performed by using two counter- 
propagating beams aligned perpendicularly to the weak 
axis of the trap. Spectra were taken by pulsing on the 
light shortly before switching off the trap and determin- 
ing the number of scattered atoms as a function of the 
frequency difference between the two Bragg beams. Since 
the kinetic energy of the scattered atoms was much larger 
than the mean-field energy, they were well separated 
from the unscattered cloud after a typical ballistic ex- 
pansion time of 20 msec. Center frequencies and widths 
were determined from Gaussian fits to the spectra. 

Duration, intensity and detuning of the Bragg pulses 
had to be chosen carefully. The instrumental resolu- 
tion is limited by the pulse duration St pu i se due to its 
Fourier spectrum, in our case requiring St pu i se > 250/xs 
for sub-kHz resolution. The maximum pulse duration 
was limited to less than one quarter of the trap period, 
Stpuise < 500^s, by which time the initially scattered 
atoms would come to rest and thus would be indistin- 
guishable from the unscattered atoms in time-of-flight. 
The light intensity was adjusted to set the peak efficiency 
to about 20 %. Sufficient detuning was necessary to avoid 
heating of the sample. The ratio of the two-photon rate 
cjfj./4A to the spontaneous scattering rate ui 2 i r/2A 2 is 
T/2A, where ujr denotes the single beam Rabi frequency, 
A the detuning and T the natural linewidth. Sponta- 
neous scattering was negligible for the chosen detuning 
of 1.77 GHz below the 3S 1/2 F=l -> 3P 3/2 F=2 transi- 
tion. 

The relative detuning of the two Bragg beams was re- 
alized in two ways. In one scheme, a beam was split 
and sent through two independent acousto-optical mod- 
ulators driven with the appropriate difference frequency, 
and then overlapped in a counterpropagating configura- 
tion. Alternatively, a single beam was modulated with 
two frequencies separated by the relative detuning and 
backreflected. Both methods are insensitive to frequency 
drifts of the laser since the Bragg process only depends 
on the relative frequency of the two beams, which is con- 
trolled by rf-synthesizers. The second method simulta- 
neously scattered atoms in the +x and —x directions and 
was thus helpful to identify motion of the cloud. We es- 
timate that residual vibrational noise broadened the line 
by less than 1 kHz. This resolution corresponds to a ve- 
locity resolution of 0.3 mm/s or 1 % of the single-photon 
recoil. At a radial trapping frequency of 200 Hz, we had 
to avoid any motion of the cloud with an amplitude larger 
than 0.2 \xm. 

Fig. 1 shows typical spectra, taken both for a trapped 
condensate and after 3 ms time of flight when the mean- 
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field energy was fully converted into kinetic energy. The 
rms-width of the resonance for the ballistically expand- 
ing cloud is 20 kHz, which is much narrower than the 
65 kHz wide distribution of a thermal cloud at 1 /iK, a 
typical value for the BEC transition temperature under 
our conditions. 




Frequency v-v [kHz] 

FIG. 1. Bragg resonances for a trapped condensate (circles) 
and after 3 ms time of flight (triangles). For comparison, the 
momentum distributions of the ground state of the trapping 
potential (dots) and of a 1 /iK cold, thermal cloud (dashes) 
are indicated. 

We could not measure the thermal distribution with 
the same pulse duration as for the condensate since the 
fraction of scattered atoms was too small due to the broad 
resonance. The spectra for the thermal cloud and the 
expanding condensate correspond to the spatial distri- 
butions observed by absorption imaging after sufficiently 
long time of flight. With this technique, the BEC transi- 
tion is indicated by a sudden narrowing of the time— of- 
flight distribution by a factor of three. Using Bragg spec- 
troscopy, the signature of BEC is much more dramatic 
- the condensate resonance is more than thirty times 
narrower than of the thermal cloud, and indeed narrower 
than the ground state of the trap. This demonstrates 
that Bragg spectroscopy is a superior way to probe the 
BEC transition. 
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FIG. 2. Mean-field acceleration of a condensate released 
from the trap. Shown is the increase of the width of the Bragg 
resonance during the expansion. The solid line is the theoret- 
ical prediction [18] using the trap frequency v T =195 Hz. The 
dashed line represents the contributions of mean-field energy 
and finite size to the total width. 

Fig. 2 shows the conversion of mean-field energy into 
kinetic energy after the atoms are released from the trap, 
indicated by a broadening of the Bragg resonance from 
about 2 kHz to 20 kHz. After 3 ms, the cloud has reached 



its asymptotic velocity Voo . Fig. 2 agrees with the veloc- 
ity evolution v r = v 00 2Tn' r t/ y/1 + (2nv r t) 2 as expected 
from the scaling laws given in ref. |L8| which are derived 
for cigar shaped condensates with large aspect ratios in 
the Thomas-Fermi approximation. In order to compare 
with the data, the Doppler width due to the velocity 
evolution, the mean-field width and the finite-size width 
were added in quadrature, assuming Gaussian shapes. 
The finite-size width was calculated from the predicted 



evolution of the size x r = xoy/l + (2irv r t) 2 |18|. 

The narrow resonance of the trapped condensate (Fig. 
1) was studied as a function of the condensate density 
and size. Fig. 3 (a) demonstrates the linear dependence 
of the frequency shift on the density. The slope of the 
linear fit corresponds to (0.54 ± 0.07) n U / h, in excellent 
agreement with the prediction of 4noU/7h. In Fig. 3 
(b), the expected widths due to the mean-field energy 
and finite size are shown for the two different trapping 
frequencies studied. The data agree well with the solid 
lines, which represent the quadrature sum of the two con- 
tributions. To demonstrate the finite-size effect the same 
data are shown in Fig. 3 (c) after subtracting the mean- 
field broadening and the finite pulse-length broadening 
(0.5 kHz). The linewidths are consistent with the ex- 
pected 1/xq dependence. Even without these corrections 
the measured linewidths are within 20 % of the value 
expected due to the Heisenberg-uncertainty limited mo- 
mentum distribution (Fig. 3 (b)). 

The momentum spread of the condensate is limited by 
its coherence length x c which, in the case of long-range 
order, should be equal to the size xq of the condensate. 
Our results show that x c f=a xq in the radial direction of 
the trap. This quantitatively confirms the earlier quali- 
tative conclusion reached by interfering two condensates 
[jl"o|l - In particular, our measurements indicate that the 
condensate does not have phase fluctuations, i.e. that it 
does not consist of smaller quasi-condensates with ran- 
dom relative phases. It would be interesting to study this 
aspect during the formation of the condensate, e.g. after 
suddenly quenching the system below the BEC transi- 
tion temperature |20| , and observe the disappearance of 
phase fluctuations. 

In this work we have determined the normal mode of 
the condensate at a momentum transfer of two photon 
recoils, corresponding to a energy transfer vq = 100 kHz. 
Different momentum transfers are possible by changing 
the angle between the Bragg beams and/or the order 
N, thus enabling measurements of the dynamical struc- 
ture function S'(q, v) over a wide range of parameters. 
At low momentum transfer, the lineshape is dominated 
by the mean-field energy and by phonon-like collective 
excitations, whereas at high momentum transfers, the 
linewidth mainly reflects the momentum distribution of 
individual atoms. This is analogous to neutron scat- 
tering in liquid helium, where slow neutrons were used 
to observe the phonon and roton part of the dispersion 
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curve, and fast neutrons were used to determine the zero- 
momentum peak of the condensate 0. While we have 
observed higher-order Bragg scattering up to third order 
in the trapped condensate using higher laser intensities, 
its spectroscopic use was precluded by severe Rayleigh 
scattering, and would require larger detuning from the 
atomic resonance. 
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FIG. 3. Bragg spectroscopy of a trapped condensate. Line 
shifts and widths are shown for various densities and sizes of 
the condensate using two different radial trapping frequencies, 
u r = (195±20) Hz (circles), and v r = (95±20) Hz (triangles). 
The lines in b) show the contributions of the mean-field en- 
ergy (dashed), due to the finite size (dotted, for both trapping 
frequencies), and their quadrature sum (solid lines). Fig c) 
displays the width after subtraction of the contributions of 
the mean-field and the finite pulse duration and compares it 
with the prediction for the momentum uncertainty due to the 
finite size. The error bars are la errors of the Gaussian fits 
to the data. 

The use of inelastic light scattering to determine the 
structure factor of a Bose-Einstein condensate was dis- 
cussed in several theoretical papers [||^] . It would require 
the analysis of scattered light with kHz resolution and 
suffers from a strong background of coherently scattered 
light H ■ Bragg spectroscopy has distinct advantages be- 
cause it is a stimulated, background-free process in which 
the momentum transfer and energy are pre-determined 
by the laser beams rather than post-determined from 
measurements of the momentum and energy of the scat- 



tered particle. 

In conclusion, we have established Bragg spectroscopy 
as a new tool to measure properties of a condensate with 
spectroscopic precision. We have demonstrated its ca- 
pability to perform high-resolution velocimetry by re- 
solving the narrow momentum distribution of a trapped 
condensate and by observing the acceleration phase in 
ballistic expansion. Since the momentum transfer can be 
adjusted over a wide range, Bragg spectroscopy can be 
used to probe such diverse properties as collective excita- 
tions, mean-field energies, coherence properties, vortices 
or persistent currents. 
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